First and second sound modes at finite temperature in trapped Fermi gases from BCS to BEC 
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We determine the temperature T dependence of first and second sound mode frequencies for trapped Fermi 
gases undergoing BCS to Bose Einstein condensation (BEC) crossover. Our results are based on the two fluid 
equations in conjunction with a microscopic calculation of thermodynamical variables. As in experiment and at 
unitarity, we show that the lowest radial breathing mode is T independent. At finite T, higher order breathing 
modes strongly mix with second sound. Their complex T dependence should provide an alternative way of 
measuring the transition temperature, T c . 
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The recent discovery of the superfiuid phases of trapped 
Fermi gases has led to considerable interest in their collective 
mode spectrum Among the modes of ex- 

perimental interest are breathing modes as well as propaga ting 
first sound. While originally theoretical attention ]8|, |9l [l0(] 
was focused on ground state properties experimental mea- 
surements are naturally not restricted to temperature T = 0. 
Indeed there is an interesting body of information which is 
emerging in these Fermi gases about the finite temperature 
behavior S 0] of the breathing modes and, more recently 
about the propagating sound velocity |2D. 

The purpose of this paper is to compute sound mode fre- 
quencies in spherically trapped Fermi gases undergoing BCS 
to Bose-Einstein condensation (BEC) crossover, at general T. 
We present a solution of the linearized two fluid equations 
and compare with recent experiments. We focus on the radial 
breathing modes and present predictions for second sound, as 
well. The structure of the two fluid equations for Bose 1 1 10 
and Fermi gases 111 211 has been rather extensively discussed. 
In the crossover regime, the normal fluid is novel lfl3i [l4ll . 
containing both fermions and noncondensed pairs, which have 
not been systematically included in previous collective mode 
literature. 

Of great importance to the field as a whole is the future 
possibility of second sound observations, in large part because 
this may ultimately help assign more precise experimental val- 
ues to the transition temperature T c . While existing experi- 
mental techniques such as vortex observation II 1 511 and den- 
sity profile features II 1 611 help establish superfluidity they pro- 
vide lower bounds on T c or determine its value for the special 
case of a population imbalanced system. Thermodynamical 
experiments measure T c more directly |fl7l[l8ll but have been 
confined to unitarity. Thus other techniques, such as second 
sound observation will be of great value. One of the princi- 
pal results of the present paper is an analysis of how second 
sound is coupled to the breathing modes. We demonstrate that 
higher order breathing modes will reveal T c through this cou- 
pling, and therefore are an alternative to direct second sound 
measurements. However, the lowest breathing mode appears 
to be remarkably T independent at unitarity. This has been ob- 
served experimentally fl, \^ and argued to follow from isen- 



tropic considerations l6fl. Here we show, that even when we 
treat the full coupling between first and second sound, we ob- 
tain similar T independent behavior at unitarity. 

At the core of the two fluid theory is the assumption that 
hydrodynamics is valid and that there are frequent collisions 
which produce a state of local thermodynamical equilibrium. 
Although there are some exceptions 11911 . reaching the two 
fluid regime has not been easy for atomic Bose gases. Two 
fluid dynamics are more readily achieved for Fermi gases, 
principally because in the crossover regime the large scatter- 
ing lengths produce sufficient collisions. Nevertheless there 
has been considerable theoretical interest in setting up 11211 
and solving B20I1 the two fluid equations for Bose conden- 
sates. Indeed, hydrodynamical approaches have successfully 
addressed both the T = and normal state regimes of the 
Bose gases [21] and demonstrated that the breathing mode 
frequencies are the same. Here, by contrast we address the 
Fermi gas case in a trap. Because they interact more strongly 
near unitarity, hydrodynamical descriptions have been argued 
quite convincingly J^MH^l to be valid. 

Previous theoretical work has been confined to T = treat- 



ments of a harmonic trap, or to finite T theories 112211 of a 
uniform gas. Our work is most similar in spirit to an earlier 
Bose gas study [20] although we introduce different numeri- 
cal techniques as well as address fermions rather than bosons. 
We note that the input thermodynamics of systems undergo- 
ing the BCS-BEC crossover which is used in the present pa- 
per, has been rather well calibrated against experimental mea- 
surements in Ref. lfl7ll and is based on a finite temperature 
extension of the simplest (BCS-Leggett) ground state. In the 
absence of a trap our results are for the most part similar to 
those in Ref. 12211 . 

We begin with the two fluid equations which describe the 
dynamical coupling of the superfiuid velocity v s and the nor- 
mal velocity v„. Just as in the spirit of the original Landau 
two fluid equations we ignore viscosity terms. In the presence 
of a trap potential V ext = \mu ho r 2 , the two fluid equations 

V(/x + Vext + m^)=0, l+v-n = 



are given by 
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j = n s v s + n„v„. Here fi is the chemical potential, p the 
pressure, and s the entropy /?<?r particle. Moreover, we have 
n n' v n + n s v s = fiv. Here n s (v s ) and n n (v n ) denote the 
superfluid and normal densities (velocities), respectively. We 
use the subscript "0" to denote equilibrium quantities such that 
v s0 = v„ = 0, V(/U + V ext ) = 0, Vp = -n VV exU and 
no, so, /io, Pa are independent of time t. Combining this with 
the thermodynamic relation du = —sdT + dp/n, we have 
VTo = 0, implying that temperature Tq is constant in the trap. 
It then follows that in equilibrium u = /io — V ex t, consistent 
with the Thomas-Fermi approximation. 

For small deviations from equilibrium, we may linearize the 
two fluid equations. Eliminating the velocities v s and v n , one 
finds Eii 
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We will focus on <5/i(r) and 5T(r) as the principal variables. 
This choice, which is different from that in Ref. 12011 . is made 
because both variables are non vanishing at the trap edge so 
that in a basis set expansion they will satisfy the same bound- 
ary conditions. Moreover, the two fluid equations are simplest 
in this form. Expressing Ss, Sp and Sn in terms of Sn, ST, the 
two fluid equations can be rewritten as 
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It is understood that all coefficients of Sn and ST are calcu- 
lated in equilibrium so that we drop the subscript "0". The 
thermodynamical quantities in equilibrium can be calculated 
following Ref. (23J, based on the standard local density ap- 
proximation (LDA), ix(r) = /i — V ext (r). Their derivatives 
with respect to T and /.i can be calculated analytically, and 

their gradients can be obtained via V/ = W el (, 

where / denotes any of the variables (n, n s ,n n , s). 

To solve the two coupled differential equations (0 and 
we assume a simple harmonic time dependence <5/i, 5T oc 
e -iwt ^ e cast jjjg differential two fluid equations into an 
eigenfunction problem with oj 2 playing the role of eigenvalue 
and the eigenfunctions given by the amplitudes of 5[i and ST. 



Since neither T nor (i depends on the density, they will not 
vanish at the trap edge. Our boundary conditions require that 
all thermodynamic variables be smooth (but not necessarily 
zero) at the trap edge. At finite T, the density in the trap de- 
creases exponentially when the local chemical potential be- 
comes negative at large radius. We choose, thus, to expand 
<5/i and ST in terms of Jacobi polynomials. For our numerics 
we choose the matrix dimension to be 300; we have similarly 
investigated matrices of dimension 200 up to 900, and found 
little change in our principal findings. 

We now turn to an important aspect of our numerics. Be- 
cause we generate some 300 frequencies in our numerical ap- 
proach it is essential to establish a mechanism for systemati- 
cally identifying first and second sound modes. To help find 
such a "fingerprint", we introduce a "decoupling approxima- 
tion" based on reducing Eqs. (01 and (|4]i to <5/i = gi(5/j,, 6T), 
5T = cj2 (S[x, ST), where gi } 2 are known functions. We elim- 
inate cross terms by setting ST = in g\ and = Q in 
<72- With these two decoupled equations it is then relatively 
straightforward to associate a profile plot of the numerically 
calculated Sp(r) and nSs(r) vs r within a trap with first or sec- 
ond sound-like modes. Here we convert S[i and ST to Sp and 
Ss via Sp = nSfi + nsST and Ss = 5fx + ST. 

We stress that this procedure differs somewhat from that in 
Ref. lf20ll where it was argued that one could associate second 
sound with a mode in which the oscillations in the thermody- 
namical variable 8T(r) were constrained to be in the conden- 
sate region of the trap. As explained below, we do not find 
this to be the case. An important check on our procedure is 
that we find that there is no sign of second sound above T c . 

Up to this point everything is general applying to both 
Fermi and Bose superfluids. All that is needed is a micro- 
scopic theory for thermodynamical variables. Here we use a 
calculational framework we have developed for treating BCS- 



BEC crossover in trapped Fermi gases H23L 12411 which em 



phasizes the importance of pseudogap effects or finite mo- 
mentum pairs. The local thermodynamical potential (density) 
f2 = Clf + Oft is associated with a contribution from gapped 
fermionic excitations ilf as well as from non-condensed pairs, 
called SI;,. We have 

n f = -— + Y / [(^-Eu)-2T\n(l + e- E ^ T )}, 

9 Vr 



-ZA 2 fi pair + 2jrin(l - e 



f! q /T 



(5) 



Here [i va i r is the chemical potential of the pairs which is 
zero below T c and the pair dispersion 57 q , along with the (in- 
verse) residue Z, can be derived from a microscopic T-matrix 
theory, described elsewhere lfl3l [l4ll . Using Q one then ar- 
rives at thermodynamical properties such as the entropy den- 
sity ns — — |p as well as self consistent equations for the 
total excitation gap A, the contribution to A from noncon- 
densed pairs (called the pseudogap) and the number equa- 
tions. These self-consistent (local) equations are simply given 
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Figure 1: (Color online) Behavior of the first (upper row) and second 
sound (lower row) modes for a spherical trap at unitarity within the 
"decoupling approximation" (see text). The left column shows the 
T-dependence of the frequencies, while the right shows correspond- 
ing typical spatial oscillation profiles for Sp [(black) solid lines] and 
nSs [(red) dashed lines] at the lowest frequencies, which provide fin- 
gerprints of first and second sound. For the second sound, nSs dom- 
inates and changes sign within the condensate. Here T c ~ 0.27Tf, 
and the arrows indicate the condensate edge. 
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number constraint N = J cPrn(r). When T < T c , we can 
use the gap equation and chain rule to eliminate the variable 

a 2 n f , d 2 n f 



A, which is a function of /i and T, via ^ = — g^g^ / 1 



e 2 n^_ / d 2 n f 



Similarly, when T > T c we use 
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the gap and pseudogap equations to eliminate the variables A 
and n P air to arrive at thermodynamical quantities. 

Figure 1 shows the lowest two collective modes at unitar- 
ity in a spherical trap, obtained by solving the chemical po- 
tential or temperature fluctuation equation in the decoupling 
approximation scheme described earlier. It is evident that our 
approximated breathing mode frequency [Fig. 1(a)] is inde- 
pendent of temperature. We understand this result by noting 
that the decoupled equation for the breathing mode is given by 
-lo 2 S(i = C MM i y 2 Sfi + C w2 • VSfi where = n(f£j 

and C PM 2 = V/i. At unitarity, n(jS^ — The only 

T-dependence contained in /irj = n{r — 0) can be elimi- 
nated via a simple re-scaling of r — > r ^/ j mus\ , yielding 
a T-independent breathing mode frequency. These arguments 
can be shown to be equivalent to the isentropic assumption of 
Ref. D6D, where v s = v„ is assumed in the two fluid model, 
leading to two simpler coupled equations associated with the 
Euler and the continuity equations. 

By contrast, the second sound mode frequency we obtain 
increases rapidly with temperature. Some typical oscillation 
profiles of Sp(r) and nSs(r) are shown in the right two pan- 
els of Figure 1 . Although the "entropy density oscillations" 
n5s(r) fall off at large r, the entropy per particle 5s(r) os- 
cillations (not shown) increase very rapidly upon entering the 
normal region. Consequently temperature fluctuations ST be- 




Figure 2: (Color online) Typical spatial oscillation profiles for Sp 
[(black) solid lines] and nSs [(red) dashed lines] obtained from the 
fully coupled equations for a spherical trap at (a-c) unitarity and (d) 
I/Ufcl — 1 at different T, for the first (top row) and second (bottom 
row) sound modes. Also shown in (a) is the T — analytical result 
[(green) dot-dashed curve]. The arrows indicate condensate edge. 



come large at the trap edge. 

Our identification of first sound for the decoupled case leads 
us to associate this mode in a coupled situation with a profile 
for which at the trap center Sp has a large amplitude, while 
nSs is almost zero (with a small peak near the trap edge) as 
in Fig. 1(b). By contrast, in the trap center the second sound 
mode has large entropy fluctuations, while the pressure fluctu- 
ations are almost zero [Fig. 1(d)]. These features will serve as 
fingerprints for distinguishing (lower order) first and second 
sound modes at finite temperatures. 

Figure 2 shows some typical eigenfunction profiles of the 
lowest modes obtained in the spherical trap upon solution of 
the fully coupled two fluid equations. The first row corre- 
sponds to the breathing mode in the unitary case. The good 
agreement between the very low T result for Sp and the T = 
analytical solution (green dashed line) in Fig. 2(a) helps vali- 
date our numerical scheme. Figure 2(b) corresponds to a high 
temperature breathing mode. In this regime the pseudogap re- 
gion outside the superfluid core is relatively large and the peak 
in n8s(r) is accordingly very broad. The lower two panels 
correspond to second sound modes for the unitary and BEC 
cases. By contrast with the breathing modes, here nSs has a 
larger amplitude than Sp with an opposite sign. Clearly this 
is very similar to what we observed in the decoupled mode 
analysis of Fig. 1. 

Figures 3(a-c) address the fully coupled equations and show 
the behavior of the lowest breathing (upper branch) and sec- 
ond sound mode (lower branch) frequencies as a function of 
temperature in a spherical trap for 1/kpa = 1, 0, and —0.5, 
respectively. For all three values of 1/k fcl we find very little 
sign of T c in the lowest breathing mode frequency. In con- 
trast, the second sound mode frequencies increase with T and 
disappear above T c . Figure 3(d) presents a more complete 
series of modes for the unitary case. Here the lines serve as 
guides to the eye for the lowest (blue) and higher order (green) 
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Figure 3: (Color online) Temperature dependence of breathing mode 
and second sound frequencies in a spherical trap. Panels (a-c) are for 
the near-BEC, unitary, and near-BCS cases, respectively. The upper 
and lower branches in (a-c) represent the lowest breathing mode and 
second sound frequencies, respectively. In (d), more complete results 
(open circles) are shown at unitarity, where the lines serve as guides 
to the eye for the breathing mode (nearly horizontal blue and green 
lines) and second sound (increasing red curves) frequencies. 



breathing mode, and second sound (red) frequencies. 

At unitarity [Fig. 3(d)] one can identify a sequence of 
higher order breathing modes which precisely overlap with 
analytical calculations for T = 0. Importantly, only the low- 
est of these is found to be a constant in temperature; the others 
are found to mix with second sound modes, as indicated in the 
solid and dashed (red) lines in Fig. 3(d). The behavior of the 
lowest mode helps to justify the isentropic assumption made 
in Ref. 0]. We understand this by referring back to the de- 
coupled profiles in Figs. 1(b) and 1(c), which are seen to be 
quite distinct. By contrast the profiles of the decoupled first 
and second sound modes at higher order (not shown) appear 
more similar to each other than their lowest order counter- 
parts. Indeed, the behavior of the profiles at higher order is 
associated with an increasing number of nodes in the curves 
of Sp and nSs (not shown), which leads to a greater similarity 
between first and second sound profiles and helps explain why 
the higher order modes are more readily mixed B25I1 . 

In summary, we have presented predictions for future ex- 
periments on higher breathing modes and second sound in a 
trap. We find that only the lowest breathing mode frequency 
has very weak T dependence. For the unitary case this tem- 
perature insensitivity was clearly observed by Thomas and co- 
workers [6]. As a result of this experiment it should not be 
surprising that we find relatively weak dependencies on ei- 
ther side of the Feshbach resonance for this breathing mode 
frequency. In the literature there are experimental claims (at 
I/Ufo, = 1.0) which are consistent with a decrease J3l in the 
radial breathing mode frequency, as indeed we find here, al- 
though ours is probably too weak to reconcile the different 
findings in Refs. [4] and J^Ht]. Finally, our more complete 
studies at unitarity show that if higher order radial breathing 
modes can be accessed, because of their strong hybridization 



with second sound, it may be possible to use these breathing 
modes, rather than direct second sound to determine the tran- 
sition temperatures T c . This should be of value since there are 
currently few experiments which can assign a value to T c over 
the wider crossover regime. 
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